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MSSM Higgs bosons are the most promising way to discover Higgs physics at hadronic col- 
liders since their cross section is enhanced compared to that of the Standard Model. We will 
present theoretical predictions for their production and decay in the Higgs— > tt channel at 
the Tevatron and the LHC, focusing on the theoretical uncertainties that affect them. The 
inferred SUSY Higgs bounds on the [tan/?; Ma] plane and the impact of these uncertainties 
will also be discussed. 



1 Introduction 

The search for the Higgs bosons which are a trace of the electroweak symmetry breaking^! is 
the main goal for current high-energy colliders. In the minimal super symmetric extension of 
the Standard Model (SM), one of the most attractive solutions of the hierarchy problem in the 
two Higgs doublets are required to cancel anomalies, which then lead to five Higgs states: 
the CP-even h, H, the CP-odd A and the two charged Higgs bosons H^. 

At tree-level two parameters in the Minimal Super symmetric Standard Model (MSSM) 
describe the Higgs sector: the vacuum expectation values (vev) ratio between the two Higgs 

doublets tan/3 = — and the CP-odd Higgs mass M4. 

At high tan j3 values, tan (3 > 10, either h or H is SM-like and its couplings to other particles 
are the same as those of the SM Higgs boson, while the other CP-even state behaves as the 
CP-odd A: same couplings and almost same mass. We will denote these two states as $ in the 
next sections. This behaviour occurs in current MSSM Higgs benchmarks scenarios^ which are 
considered at the Fermilab TevatrorP and the CERN LH(P3 colliders. 

6-processes are dominant as they are proportionnal to tan f3 contrary to that of the top-lo op. 
We will thus consider t he gluon-gl uon fusion Higgs production through bottom quark loorP^I 
and the 66 fusion channeff^HECniS followed by the Higgs— > t + t~ desintegration. Squark loops 
can be safely neglected while SUSY corrections to the $66 coupling nearly cancel out in the 



production cross section times branching ratio calculation^ 

We will present numerical results at the Tevatron and the 1HC (LHC at 7 TeV) for tan f3 = 1, 
which means that we have to multiply by 2 tan 2 /3 for actual values. Theoretical uncertainties will 
also be presented and their implications on the MSSM parameter space limits will be discussed. 
A more detailed discussion can be found in 
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Figure 1: o-g^j;^ and cr£r_^$ central cross sections using MSTW 2008 PDFs and &bb unit couplings together with 
total uncertainties at the Tevatron (left) and at the 1HC (right). In the insert are shown the individual sources of 

uncertainties normalised to the central cross section. 



2 SUSY Neutral Higgs production at the Tevatron and the 1HC 



2.1 gg — > $ channel 

The Higgs bosons in the gluon-gluon fusion channel is produced through top and bottom quarks 
loops. At tan (3 > 10 values the top loop is strongly suppressed because &ti is inversely propor- 
tionnal to tan /3 contrary to the $66 coupling. Although the top loop is known up to next-to- 
next-to-leading order (NNLO) in QCD, the 6-loop is known up to next -to-leading order (NLO) 
onlJS. We will use NLO MSTW 2008 parton distribution functions (PDF) sePwe consider 
the standard QCD theoretical uncertainties that have been discussed in 

It is customary to estimate the uncertainty due to the missing higher order terms in a 
perturbative calculation by varying the renormalisation scale [jlr and the factorisation scale \xp 
around a central scale /ir> ^ < Hr,^f < ^Mo- We take /io = ^M<j, as the central scale in order 
to be consistent with the SM calculatiorff^, as one of the CP-even Higgs is SM-like. We use 
K = 2 in the gluon-gluon fusion channel and obtain Aa/a ~ ±20% at the Tevatron (±15% at 
the 1HC). 

The next source of uncertainties is due to the j^ombined uncertainty from the PDF and 
as coupling. We use MSTW collaboration scheme^ and calculate the PDF±A exp+th as 90% 
CL uncertainty which is equivalent to the MSTW PDF4LHC recommandatiorP^ an d we obtain 
Aa/a ~ ±10% both at the Tevatron and the 1HC. 

The last important uncertainty is specific to the MSSM case, and deals with the mj, mass. 
There are two types of uncertainties: the experimental errors on the MS m^m^) value and the 
uncertainty due to the scheme choice for the renormalisation of the 6-mass. The first uncertainty 
will cancel out in the production cross section times branching ratio (see below) but not the 
other one. We obtain Aa/a ~ ±15% at both colliders due to these 6-quark issues. 

All these individual sources of uncertainties are shown in Fig{TJ in the insert. We also 
display the total uncertainty on the cross section when combining the uncertainties according 
to the procedure developed in Refill We obtain Aa/a ~ ±58%,— 40% at the Tevatron and 
Aa/a ~ ±53%, -38% at the 1HC. 



2.2 bb -> $ channel 

The bottom quark fusion channel is strongly enhanced because of the tan f3 effect in 6-quark 
processes. This channel is known in the SM up to NNLO in QClP^and we rescale the predictions 
with the MSSM $66 coupling to obtain a NNLO MSSM prediction. We use the same PDF set 




Figure 2: a(pp — > A) x BR(y4 — > t + t~) as a function of Ma at the Tevatron (left) and at the 1HC (right), together 
with the associated overall theoretical uncertainty; the uncertainty when excluding that on the branching ratio is 
also displayed. In the inserts, shown are the relative deviations from the central values. 



as for gluon-gluon fusion and consider the same set of theoretical uncertainties. 

For the scale uncertainties we consider here k = 3 instead of k = 2 in gg — > <£. Indeed this 
is well known that either four or five active flavours schemes can be used for the calculation. 
The two predictions differ significantly^ and one way to reconcile them is to allow such a scale 
interval. Furthermore this also allows the inclusion of the 6-mass scheme uncertainty that was 
obtained separately in the gluon-gluon fusion calculation. We obtain in the end Aa/a ~ 30% 
for low masses at the Tevatron (±25% at the 1HC). 

The combined PDF+ces uncertainty is calculated exactly as in the gluon-gluon fusion case. 
We obtain in the bottom quark fusion Aa/a ~ ±20% for low masses and ~ ±30% for high 
masses at the Tevatron (~ 10% at low masses and ~ ±20% at high masses at the 1HC). 

The last uncertainty involves only the experimental 6-mass error. We obtain a ±10%, —4% 
uncertainty at the Tevatron (nearly the same at the 1HC), which as discussed in the next section 
will cancel out in the cross section times branching ratio calculation. 

All the uncertainties are displayed in FigfTJ The overall total uncertainty is Aa/a ~ 
±50%, -40% at the Tevatron (±40%, -30% at the 1HC). 



2.3 Combinaison with the <3? — > tt branching ratio 

We finally evaluate the combinaison of the two production channels together with the branching 

rati o <E ' — » t + t~. The issue is how to combine the uncertainties and we proceed as stated in 
RefsPEl 

the cross section uncertainties are weighted according to their importance and we 
add linearly the decay branching ratio uncertainty which is ~ ±4%, —9% on BR($ — > t + t~ ~ 
10%jP^. In this procedure, as the uncertainties due to the experimental errors on 6-mass are 
anti-correlated in the production and decay, they cancel out. 

We then obtain A(a x BR)/(<r x BR) ~ ±50%, -39% at the Tevatron and ~ ±35%, -30% 
at the 1HC, as shown in Figj2j 



3 Higgs bounds on the MSSM parameter space 

We are left to evaluate the impact of the theoretical uncertainties calculated above on the 95% 
CL limits in the [tan /?; Ma] plane using the experimental results at the Tevatron and the 1HC. 
The results presented above are quite model independant as they do not depend on the details of 
the MSSM model as long as we have a degeneracy in the h/H,A spectrum. We apply the limit 
on the minimal cross section times branching ratio instead of the central prediction in order to 
take into account the theoretical uncertainties. 

The result is shown in Figl3]and the theoretical uncertainties are extremely important. We 
obtain tan j3 > 45 at the Tevatron, which thus reopens a large part of the parameter space 




Figure 3: Contours for the expected a(pp — > $ — > t + t rate at the Tevatron (left) and at the 1HC (right) in the 
[Ma; tan /3] plane with the associated theory uncertainties, confronted to the 95% CL exclusion limit obtained by 

CDF/DO and CMS. 



excluded by CDF/DCPJ The comparaison with CMS results at the 1HC shows a slight reduction 
of the exclusion limit as we obtain tan j3 > 29 to be compared with tan /? > 23. The result is 
comparable to what can be obtained with the theory uncertainty quoted by 
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